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We study an RF SQUID, in which a section of the loop is a freely suspended beam that is allowed 
to oscillate mechanically. The coupling between the RF SQUID and the mechanical resonator 
originates from the dependence of the total magnetic flux threading the loop on the displacement 
of the resonator. Motion of the latter affects the visibility of Rabi oscillations between the two 
lowest energy states of the RF SQUID. We address the feasibility of experimental observation of 
decoherence and recoherence, namely decay and rise of the visibility, in such a system. 

PACS numbers: 03.65.Yz, 85.25.Dq 



I. INTRODUCTION 

Decoherence occurs when a quantum system is cou- 
pled to a noisy environment at a finite temperature. De- 
coherence is commonly quantified by a visibility factor, 
which characterizes the relative amplitude of a measured 
interference signal. In many cases the main contribu- 
tion to decoherence originates from the many degrees 
of freedom of the environment, which all have a simi- 
lar coupling strength to the interfering degree of freedom 
of the quantum system. In such a case the visibility fac- 
tor is expected to decay monotonically as a function of 
time (typically, the decay is exponential). On the other 
hand, when only a few degrees of freedom in the envi- 
ronment significantly contribute, the time dependence of 
the visibility factor is not necessarily monotonic. Reco- 
herence occurs when the visibility factor increases with 
time. Experimental demonstration of this phenomenon 
is important since it may provide a crucial test to the 
theory of quantum measurement Q, • Decoherence and 
recoherence were recently discussed theoretically in Refs. 
H II H II II H The interfering quantum system 
in Ref. [3| was a single level quantum dot, in Refs. 
0,0,0,0 it was an optical mode in a cavity, and in Refs. 
[a, |j| a superconducting charge (Cooper-pair box) and 
phase Josephson qubit, respectively. In all these cases, 
the interfering quantum system is coupled to a vibrat- 
ing mode of a mechanical resonator (typically the lowest, 
fundamental mode). Recoherence can occur in such sys- 
tems provided that the coupling between the interfering 
quantum system and the mode of the mechanical res- 
onator is made sufficiently strong, whereas the coupling 
to other degrees of freedom in the environment is suf- 
ficiently weak. Satisfying this condition experimentally 
when the interfering degree of freedom is a single elec- 
tron, as in the Ref. or a single photon, as in Refs. 
0,1a, la, 0, turns out to be very difficult. 

In the present paper we study an alternative configu- 
ration consisting of an RF superconducting quantum in- 
terference device (SQUID) integrated with a mechanical 
resonator in the shape of a doubly clamped beam. The 



dependence of the total magnetic flux threading the loop 
on the beam's displacement leads to a coupling between 
the RF SQUID and the mechanical resonator. We study 
the effect of such a coupling on the visibility of Rabi os- 
cillations between the two lowest energy states of the RF 
SQUID, and discuss the required conditions for experi- 
mental observation of decoherence and recoherence orig- 
inating from the coupling to the mechanical resonator. 

The paper is organized as follows. The Hamiltonian 
for the closed system is obtained in section II. An adi- 
abatic approximation is employed in section III to sim- 
plify the equations of motion of the system by consider- 
ing the mechanical motion as slow in comparison with the 
faster dynamics of the RF SQUID. Further simplification 
is achieved in section IV by taking into account only the 
two lowest energy levels of the RF SQUID. In section 
V we calculate the effect of the mechanical resonator on 
the visibility of Rabi oscillations between these two en- 
ergy levels. Corrections due to finite temperature and 
mechanical damping are considered in sections VI and 
VII respectively. The validity of the adiabatic approxi- 
mation is examined in section VIII. A numerical example 
is given in section IX and discussion and conclusions are 
given in section X. 

Similar systems consisting of a SQUID integrated with 
a nanomechanical resonator have been recently studied 
theoretically. Zhou and Mizel have shown that nonlinear 
coupling between a DC SQUID and a mechanical res- 
onator can be employed for producing squeezed states of 
the mechanical resonator [lOj. More recently, Xue et al. 
have shown that a flux qubit integrated with a nanome- 
chanical resonator can form a cavity quantum electrody- 
namics system in the strong coupling region 11]. 

II. HAMILTONIAN OF THE CLOSED SYSTEM 

Consider the RF SQUID shown in the inset of Fig. ^ 
in which a section of the loop is freely suspended and 
allowed to oscillate mechanically. We assume the case 
where the fundamental mechanical mode vibrates in the 



2 




n 



1400 



1200 



1000 d 
800 O 



600 



U400 
r 



FIG. 1: (Color online) The potential U for the case 

$e = $o/2 and /3l = 20. The inset schematically shows the 
device. 



plane of the loop and denote the amplitude of this fiexural 
mode as x. Let m be the effective mass of the fundamen- 
tal mode, and loq its angular frequency. A magnetic field 
is applied perpendicularly to the plane of the loop. Let 
$ e be the externally applied flux for the case x = 0, and 
B is the component of the magnetic field normal to the 
plane of the loop at the location of the doubly clamped 
beam (it is assumed that B is constant in the region 
where the beam oscillates). The total magnetic flux $ 
threading the loop is given by 



$ = <J>„ + Blx + LI 



(1) 



where L is the self inductance of the loop, and / is an 
effective length of the beam. The contribution of other 
mechanical modes of the beam to $ is assumed to be 
negligibly small. 

A Josephson junction (JJ) having a critical current I c 
and capacitance C is integrated into the loop. We first 
consider the dynamics of the closed system consisting 
of the RF SQUID with the integrated doubly clamped 
beam. The effect of damping due to coupling to other 
degrees of freedom in the environment will be discussed 
later. 



A. Lagrangian 

The Lagrangian of the closed system is a function of 
the position x, flux $ and their time derivatives (denoted 



by overdot): 

1 C$ 2 

C = -mx 2 H U (x, $) , 

2 2 

where the potential energy U is given by 

_ mu 2 x 2 t ($ - $ e - Blx) 2 $ o/ c cos (^) 



(2) 



, (3) 



2 2L 2tt 

and <&o = h/2e is the flux quantum (see Fig. QJ. The 
resulting Euler - Lagrange equations are 



Bl 

mx + moj 2 x ($ — $ e — Blx) = , (4) 

L 



r . T , *-&e-Blx , T . /2tt$ 
C$ H ; h I c sin I -— - 

$0 







(5) 



Note that the gauge invariant phase across the Josephson 
junction 7,/ is given by 



7,7 = 27m 



2tt$ 



(6) 



where n is integer. By using this and Eq. JQ) the equa- 
tions of motion can be rewritten as 



mx + mu> x — BlI = , 



$0 

7 c sin7j + C— 7j = I 

Z7T 



(7) 



(8) 



The interpretation of these equations is straightfor- 
ward. Eq. (JJ expresses Newton's 2nd law where the 
force is composed of the restoring mechanical force and 
the Lorentz force acting on the movable beam. Whereas 
Eq. (JSJ) states that the circulating current / equals the 
sum of the current I c sin 7j through the J J and the cur- 
rent CV through the capacitor, where the voltage V is 
given by the second Josephson equation V = ($o/27r) jj. 



B. Hamiltonian 

The variables canonically conjugate to x and $ are 
p = dC/dx = mx and Q = d£/d& — C<i> respectively. 
The Hamiltonian is given by 



n 



p 

2m 



91 

2C 



U(x,$) 



(9) 

Quantization is achieved by regarding the variables x, 
p, <& and Q as Hermitian operators satisfying the fol- 
lowing commutation relations [x,p] — [&,Q] = ih and 
[ar,*] = [ar J Q] = [ Pj *] = [p,Q] = 0. 
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III. ADIABATIC CASE 

The Hamiltonian JO} can be written as H = Ho + Hi , 
where 



n 2 1 

-11 P i 1 2 2 

H ° = 2m + 2™^ 



(10) 



(11) 



Using the notation Uq = <&o/87r 2 L, = Q /Bl, and 
/?i = 2ttLI c /^> , the term it can be written as 



u = Un 



— - 2/3i cos — — 



As a basis for expanding the general state of the system 
we use the solutions of the following Schrodinger equation 



Hi \n (x)) = e n (x) \n(x)) 



(13) 



where x is treated here as a parameter (rather than a 
degree of freedom). The local eigenvectors are assumed 
to be orthonormal 



(m (x) \n(x)) = 6„ 



(14) 



The wavefunctions associated with the local eigen- 
statcs 



are the solutions of the Schrodinger equation 



h 2 d 2 



+ u($;x) 



2C d<s> 2 

The total wave function is expanded as 



i> (x,$,t) = ^2t n (x,t)\n(x)) 



(15) 



(16) 



(17) 



In the adiabatic approximation [12\ the time evolution 
of the coefficients is governed by the following set of 
decoupled equations of motion 



2m 



+ V m (x) 



(18) 



where the adiabatic potentials V rn (x) are given by 



The validity of the adiabatic approximation will be dis- 
cussed in section VIII below. 

To numerically evaluate the eigenvalues e m {x), it 
is convenient to introduce the dimensionless variables 

27T$/$o = TT + (j>, 27T$ e /$ = 7T + 4>e, ^Tix/x^ = <j) x . 

Using this notation the Schrodinger equation (|16|) can be 
rewritten as 



(20) 



Op 2 Uq' 



where f} c = 2e 2 /CU Q ,X n , x = e n (x) /Uq, and 



U 



(f>e-4>x) +2/3lCOS</> 



(21) 



IV. TWO-LEVEL APPROXIMATION 



Consider the case where \<j> e \ <C 1 (namely, $ e ~ 
$o/2), \4> x \ <C 1) and 0l > 1- In this case the local 
potential u (</>) given by Eq. 1)2 111 contains two wells sep- 
arated by a barrier near = 0. At low temperatures only 
the two lowest energy levels contribute. In this limit the 
local Hamiltonian Hi can be expressed in the basis of the 
states \~r\) and \rx), representing localized states in the 
left and right well respectively having opposite circulat- 
ing currents. In this basis, Hi is represented by the 2x2 
matrix 



Hi 



V (4>e + 4>x) 

A 



A 



(22) 



The real parameters r\ and A can be determined by 
solving numerically the Schrodinger equation (|20|l . 
Using the notation 



tan 



Hi can be rewritten as 



7] (4> e + 4>x) 



(23) 



The eigenvectors and eigenenergies are denoted as 



Hl\±)=E ± \±) , 



where 



sm ; 



H; I-) 



(25) 



(26) 



V m (x) = -mOI I + E rn (x) 



(19) 



cf + A 2 . 



(27) 
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V. RABI OSCILLATIONS 

Consider the following experimental protocol for de- 
tecting Rabi oscillations between the two lowest energy 
states of the RF SQUID. The first stage consists of state 
preparation performed by applying a large constant ex- 
ternal flux <j) e . At time t = the external flux is switched 
off and the system starts oscillating. At a later time t > 
the final state of the RF SQUID is measured. 



Thus, both mechanical states associated with the RF 
SQUID states |+) and |— ) will at t = start oscillat- 
ing with different frequencies u>q\/1 + ( and ~ C 
respectively around the point x — 0. Consider the case 
where ( < 1. Using Eq. I|31l) one finds that the approx- 
imation 



|±(x)) = |±(x = 0)) = 



(34) 



A. State Preparation 

The system is first prepared in an initial state by ap- 
plying an external bias flux 4> e such that ip e 3> A/ry. In 
this limit one finds approximately |+) = |— ) = |rv), 
and s± — ±7] {4> e + 4> x ). Thus, the adiabatic potentials 
Eq. (|I9I) are given by 



V± (x) = -mwg (x ± x ) - T) 



7TXo 

Xj, 



where 



■ro = 



2nr\ 



(28) 



(29) 



Assume also the case where the temperature T is rel- 
atively low k B T < A. In this limit the RF SQUID is 
expected to occupy its ground state |^) in thermal equi- 
librium. The mechanical resonator is expected to be in a 
thermal state of the potential well V- (x) centered at xq 
[seeEq. (J2HJ]- 



B. Switching off the External Flux 

At time t = 0, the external flux <j) e is suddenly switched 
to a new value <fi e = 0. Using the notation 



c = 



7] 27TXQ 



A a;^ 

one finds to lowest order in 



(30) 



can be employed in the region 
chanical resonator oscillates. 



< xq where the me- 



C. Measuring the RF SQUID Final State 

Consider the case where the mechanical system was at 
time t = in a given state, denoted as |£o) e i with a wave 
function £o {%)■ We first calculate the time evolution for 
a given state, and later perform a thermal averaging over 
initial states. The state of the system at t = can be 
expressed as 



t[> (t = 0) = £+ (t = 0) |+ (x = 0))+£- (t = 0) |- (x - 0)> , 

(35) 

where 



U (* = 0) = ±^ (x) 



(36) 



In the last step the state of the RF SQUID is mea- 
sured. What is the probability to find the RF SQUID in 
a given state \x) at time t ? To calculate this probabil- 
ity Pi x ) (t) one has to trace out the mechanical degree of 
freedom. By using Eq. I|17|l and employing the two-level 
approximation one finds in general 



P lx) (t) = J dx \C+ Or, t) <x| + (x)) + £_ t) ( X \ - (x))\ 2 

(37) 

As an example, consider the case where \x) = \^)- 
Using Eq. I)34|l one finds 




l±) 



(31) 



(32) 



and the adiabatic potentials (|f 9fl for this case are given 
by 



V± (x) = \miol (1 ± C) x 2 ± A 



(33) 



P\r.) (<) = -+ Re y dxi + {x, t) C (x, t) . (38) 

Alternatively, using Eqs. (|I8|I and (|33|l this can be ex- 
pressed as 



P\r.) (t) = l + J Re 



u (t) exp 



2iAt 



(39) 
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where 



and 



. iH(-Qt\ ( iH(()t\ . 
vo (t) = e (&| exp L^- exp ^- |£ , 



and 



(40) 



H(0 



1 



+ o^o (! + « 



(41) 



2m ' 2' 

Eq. (|39|) indicates that the visibility of Rabi oscillations 
(occurring at angular frequency 2A/h) is diminished by 
the factor \v (t)\ (note that in general \i/q < 1). 
The Hamiltonian H$ can be written as 



h (o = yr+c^ (o) + v ( 



(42) 



where 



on = -r-'o 



TTiWo 

2ft 



(47) 



For a general c-number a, the operator I? (a) trans- 
forms the vacuum state |0) into a coherent state \a), i.e., 
D (a) |0) = | a). Using this notation one finds 



v (t) = ( D f (a ) exp - 



iC-ff (0)i 



D(ao) 



(48) 



where the brackets () represent thermal averaging. It is 
convenient to employ the coherent states diagonal repre- 
sentation (P representation) 13] of the density operator 
at thermal equilibrium 



P 

2m 



(i-v/ITc) + ^(i-v/ITc + c) 



(43) 

The Hamiltonian \/T+ (0) is associated with a har- 
monic oscillator having mass m/yl + C and a resonance 
frequency woVT+X- Assuming that £ -C 1 one can em- 
ploy the approximation 



P = 



J J d 2 aP (a) | a) 



(49) 



where d 2 a denotes infinitesimal area in the a complex 
plane, namely d 2 a = d {Re a} d {Ima}, the probability 
density P (a) is given by 



H (C) y/l + ^H (0) 



(44) 



This approximation greatly simplifies the analysis since 
annihilation and creation operators associated with both 
Hamiltonians H (£) and (— C) are common. Note 
that the time evolution generated by both Hamiltonians, 
H (C) and y/1 + (H (0), is periodic in time with the same 
period 2tt /loo^/I + £. Thus, the error introduced by this 
approximation is small even for times much longer than 
the period time, provided that the condition ( < 1 is 
satisfied. Using this approximation and keeping terms 
up to first order in £ yield 



Mt)= e(£o|exp(-^Hf) | V . 



(45) 



VI. THERMAL AVERAGING 



At finite temperature T the term vo (t) has to be cal- 
culated by averaging over a thermal distribution of initial 
states |£o) e - At times t < the mechanical resonator is 
expected to be in a thermal state of the potential well 
U_ (x) centered at Xq [Eq. (|28|l ]. It is convenient to 
express this thermal distribution using a displacement 
operator D (ao), where 



D (a) = exp 



rriLUo 
2h 



(a — a*) x — i 



1 



2hmujQ 



(a + a*)p 



P(a) 




7r (n) 



and 



1 



e huj /k B T _ I 



is the thermal occupation number. 
Thus 



v (t) = Tr 



pD 1 * (a ) exp ( - 



i(H (0)t 



= J J d 2 aP (a) 
x (a | D' (ao) e x P 

Using the identity 



iCH(0)t 



(50) 



(51) 



D(ao) 



D(a )\a) 



(52) 



D (ao) | a) = exp 



aoa 



a n a 



\a + a) 



(53) 



(46) and noting that ao is real yield 



6 



VII. EFFECT OF MECHANICAL DAMPING 



v (*) = 



—i^uot/2 
7T (n) 
—i£uot/2 



■ exp 



d 2 a exp — 



|«| 2 + £,\ao + a\ 2 



dx exp 



-(i + u + 



/oo 
exp[-(l + C)2; 2 ] , 
-oo 



*b (*) 



-iCu t/2 



exp [-^(l-e-^*)] 



and the visibility factor in this limit is given by 



M*)r 



exp 



-4a 2 , sin 2 



V 2 



(54) 



where £ = (l — e ^ Wot ) (n) and ao = cco/V (n). 

In the limit of zero temperature where (n) — » one 
finds 



(55) 



(56) 



Another case of interest is the limit of short times. The 
term \uq (t)\ is calculated to lowest order in t using Eq. 
(|48|l and perturbation theory 



M*)r 



(57) 



where 



Consider in general a mechanical resonator in a su- 
perposition of two coherent states \a±) and |aa). Cou- 
pling between the resonator and a thermal bath at tem- 
perature T induces decoherence with a rate jd given by 



7^-g-loi-oal coth— , (62) 

where luq and Q are the resonance frequency and quality 
factor respectively. 

Damping is thus expected to further diminish the vis- 
ibility of Rabi oscillations. The factor v it) is written 
as 



v (t) = u (t) v A it) 



(63) 



where Vd (*) represents the contribution of damping. 

To provide a rough estimate of the factor Vd (t) in the 
present case the c-numbers a\ and 012 are substituted by 
the thermal average values of the distributions associated 
with the |+) and | — ) states respectively Q, and thus we 
take 



ai (t) = ao exp 
«2 (t) = ao exp 
We further require that 

dv d 
dt 

and obtain 



-' 'I 1 + 1 I wo* 



-i I 1 - I ) u t 



-IdVd 



(64a) 
(64b) 

(65) 



V H = <£>t (a ) H 2 (0) D (a ))-(D^ (a Q ) H (0) D (a )) . 

(58) 

Using Eqs. and one finds 



V H = hW {a 2 Q + („)) 



(59) 



The result can be expressed in terms of a decoherence 
rate 7 m 



v d (*) = exp 



2k B T 



sin (C^o*) 
(uj t 



(66) 



Recall that the recoherence peaks, where \vq (t n )\ = 1, 
occur at times t n — 2im/C t ujQ, where n is integer [see 
Eq. (|54|l ] . Recoherence can be detected only if jd is 
sufficiently small. For the first recoherence peak at time 
t%, we have 



where 



Mi)| 2 = l-(7m*) 2 , 



1/2 



(60) 



(61) 



Vd (*i 



exp 



87ra§ hhj 

coth 

CQ 2fc B T 



(67) 



whereas for the other recoherence peaks the following 
holds 



Vd {t n ) = [V d (ti)] n 



(68) 
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In the case Hloq <§C ksT one has 



ing parameters 



Vd (ti) = exp 



4tt / x 



(69) 



where At is the thermal length 



At = 



(70) 



VIII. ADIABATIC CONDITION 

We now return to the adiabatic approximation and ex- 
amine its validity. In the adiabatic limit the off-diagonal 
terms in the set of coupled equations for the amplitudes 
^ n are considered negligibly small, and consequently no 
Zener transitions between adiabatic states occur. This 
approximation yields the set of decoupled equations (|18|l . 
To calculate the Zener transition probability to lowest 
order we consider the off diagonal elements as a pertur- 
bation. 

Consider mechanical oscillations with an amplitude Xq 
and assume the case where <f> e — 0. A Zener transition 
is most likely to occur near the times when the mechan- 
ical resonator crosses the point x = 0, namely, when the 
mechanical velocity peaks and the energy gap £+ — e_ ob- 
tains its smallest value. The probability pz that a Zener 
transition will occur per such a crossing can be calculated 
using Eq. (C25) of Ref. 



pz = exp 



Xcf, 



2x (3h rjUo 



(71) 



where (3h = hujo/Uo. 
valid when < 1. 



The adiabatic approximation is 



IX. ESTIMATION OF PARAMETERS 



Satisfying all the above mentioned conditions required 
for experimental observation of decoherence and recoher- 
ence is quite challenging. However, a careful design to- 
gether with state of the art fabrication and cryogenics 
techniques may allow experimental implementation. We 
examine below an example of a device having the follow- 



L 

C 

Ic 

m 
Bl 

Q 

T 



6.5 x 10~ n H , 
7.4 x 1(T 17 F , 
10 a« A , 
10- 16 kg, 
640 MHz , 
T x A*m , 
10 4 , 
0.05 K . 



(72a) 
(72b) 
(72c) 
(72d) 
(72e) 
(72f) 

(72g) 
(72h) 



These parameters for both the RF SQUID [lj| and for 
the nanomechanical resonator |2(j are within reach with 
present day technology. 

The chosen value of L corresponds to a circular loop 
with a radius of about 10 Aim and a wire having a cross 
section of about (0.1 A* m ) 2 > whereas the values of C and 
I c correspond to a junction having a plasma frequency of 
about 8 THz. The parameter Bl plays a crucial role in de- 
termining the coupling strength between the mechanical 
resonator and the RF SQUID. Enhancing the coupling 
can be achieved by increasing the applied magnetic field 
at the location of the mechanical resonator B. However, 
B should not exceed the superconducting critical field. 
Moreover, the externally applied magnetic field at the 
location of the J J must be kept at a much lower value in 
order to minimize an undesirable reduction in I c . This 
can be achieved by employing an appropriate design in 
which the applied field is strongly nonuniform. 

Using these values one finds 



(3c 

Ph 
Uo 
k B 



At 



1.9 , 


(73a) 


0.78 , 


(73b) 


4.8 x 10~ 4 , 


(73c) 


64K , 


(73d) 


2.1 nm , 


(73e) 


9.1 x 10 4 , 


(73f) 


9.0 x 10~ 6 nm . 


(73g) 



The values of (3l and (ic are employed for calculating 
numerically the eigenstates of Eq. (|2"U1) . Fig. [21 (a)- 
(c) shows the first 3 lowest energy states for the case 
4>e + 4>x = 0, whereas panel (d) shows the dependence of 
the energy of the two lowest energy states on <p e + <f> x . 

From these results one finds for the values of the 77 and 
A parameters in the two-level approximation to Hamil- 
tonian Hi [Eq. 



V = 2.5f/ , 
A = O.UUq 



(74a) 
(74b) 
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FIG. 2: (Color online) Eigenstates of Hi. (a)-(c) The first 3 
lowest energy states for the case (f) £ + 4> x — 0. (d) The energy 
of the two lowest states vs. d> e + d> x . 



Using these values yields 



x 


= 4.1 x l(T 6 nm , 


(75a) 


c 


= 2.5 x 1CT 4 , 


(75b) 


x 

x<$> 


= 2.o x i(r 6 , 


(75c) 




= 0.18 , 


(75d) 


(n) 


= 1.2, 


(75e) 


(Cao^o)" 1 


= 5.6 /its , 


(75f) 


4tt fx V 

CQ \x T ) 


= 1.0, 


(75g) 


2tt 


= 6.3 /xs , 


(75h) 


x$ A 2 


= 3.0 x 10 6 . 


(75i) 


Ix^fih r/Uo 



Eqs. (|75fl) and )75gfl indicate that observation of both 
decoherence and recoherence, for the case of the present 
example, is feasible, provided that the decoherence time 
of the RF SQUID due to other mechanisms is sufficiently 
long, i.e., on the order of microseconds [2l|. Moreover, 



Eq. (|75i|l ensures the validity of the adiabatic approxi- 
mation. 



X. DISCUSSION AND CONCLUSIONS 

A possible, alternative protocol to the presently con- 
sidered one for observing decoherence/recoherence phe- 
nomena is the so-called Ramsey interference experiment 
that proceeds as follows 8] : (i) At time t < 0, the state 
is prepared in the ground state |rv), identically to the 
above considered protocol by applying an external bias 
flux </> e such that <fi e ^> A/77 i ( n ) At time t = 0, the 
external flux <p e is suddenly switched to the new value 
ef> e = 0, again just as in the above protocol, but then 
after one-quarter of a Rabi oscillation period, <j) e is sud- 
denly switched back up to the same non-zero value as 
was applied during first, preparation stage; (iii) The flux 
qubit and mechanical oscillator are then left to interact 
for a certain duration with <j> e kept constant; (iv) Stage 
(ii) is repeated again; (v) The state of the qubit is read 
out. 

The effect of stage (ii) is to prepare the flux qubit in 
a state which is an equal magnitude superposition of the 
circulating current states |rv) and |^). Each of these 
states is associated with the different spatially-shifted po- 
tentials V±(x) [Eq. J2EJ)], so that during the interaction 
stage (iii) an entangled state develops between the oscil- 
lator and flux qubit, giving rise to decoherence of the re- 
duced qubit state. After one full mechanical period, the 
entanglement is undone, resulting in recoherence. The 
second, quarter Rabi period pulse enables one to probe 
the decoherence/recoherence, simply by measuring the 
probability to be in one of the measurement basis states, 
e.g., the ground state \r\). By repeating the Ramsey 
protocol many times, allowing the interaction duration 
to range over several mechanical periods, oscillations in 
the visibility are observed providing a signature of deco- 
herence / recoherence. 

The Ramsey protocol has the obvious advantage over 
the above considered protocol (where one always remains 
at the 4> e = degeneracy point during t > 0) that the 
decoherence/recoherence times are shorter by the factor 
of However, the disadvantage with the Ramsey pro- 
tocol is that qubit decoherence times are considerably re- 
duced away from the degeneracy point. The origin of the 
reduction in these two competing timescales is of course 
the same: the mechanical oscillator and flux noise cou- 
ple more strongly (i.e., linear coupling) to the circulating 
current basis states |rv) and |^) than to the eigenstate 
basis states at the degeneracy point (i.e., quadratic cou- 
pling). Depending on how the qubit decoherence rate 
varies with the externally applied flux, it may be that 
operating a small distance from the degeneracy point 
is more favorable for observing recoherence effects pi) . 
However, the resulting coupled quantum dynamics is not 
as simple to describe as at the special limiting bias points 
where the Hamiltonian Hi [Eq. Q22))] is either (approxi- 
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mately) purely diagonal or off-diagonal. 

In the present paper we have considered a flux qubit 
in the form of an RF SQUID, a system that is relatively 
simple to analyze. However, a double well potential can 
be formed only when the inductance L is sufficiently large 
and the condition /3l > 1 is satisfied. In this limit, the 
loop is relatively large and consequently large pickup of 
external flux noise results in a relatively short flux qubit 
decoherence time On the other hand, this prob- 

lem can be partly solved by employing the configuration 
of a loop having three JJs [2^. where a portion of the 
necessary total SQUID inductance is provided by the ef- 



fective inductance of the additional JJs; the three JJ su- 
perconducting loop would likely be the preferred choice 
for experimental implementation. 
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